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Title: Method and apparatus for obtaining information about the size 
distribution of macroscopic particles in a liquid 

The invention relates to a method and apparatus in which, using 
reflection of ultrasonic pulses by macroscopic particles such as oil droplets 
or sand grains in water, information about the size distribution of the 
particles in a liquid is determined. 
5 European patent application No. 0801305 describes an apparatus for 

characterizing a suspension. The apparatus generates ultrasonic pulses in a 
beam in a liquid. Reflection of a pulse on an individual particle in the liquid 
results in an echo. If the particle is at a given position, the amplitude of the 
reflection is coupled one-to-one to the size of the particle, Accordingly, in 

10 measurements on reflections of consecutive pulses on particles of different 
sizes, different amplitudes are measured. 

The known apparatus makes a histogram of counts of the number of 
times that reflections of different amplitudes occur. The histogram contains 
information about the particle size distribution, that is A the concentration of 

15 particles as a function of the size of the particles. With information about 
the particle size distribution, in turn, for instance the total concentration of 
particles in the liquid can be determined. 

European patent application No. 0801305 describes a technique to 
extract information about the size from the distribution of the amplitudes. 

20 An expression is used for the probability that reflections of different 

amplitudes occur. This expression relates the probability to the particle size 
distribution. If the reflecting particles were always at the same place in the 
beam, the number of reflections having a particular amplitude would be 
proportional to the fraction of particles having the size leading to that 

25 amplitude. 

However, the information about the particle size distribution is 
smeared across the amplitudes becHuse* felig^eflection amplitude of a 
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particle, besides being dependent on the size of the particle, also depends on 
the position of the particle in the beam. The farther the particle is off the 
center of the beam at reflection, the smaller the amplitude. The expression 
used for the probability of a reflection measurement with a particular 
5 amplitude smears the particle size distribution across the amplitude 
distribution to give expression to this effect. 

Using a ma xim um likelihood technique, parameters of the particle size 
distribution are estimated that maximize the thus expressed probability of 
the actually measured numbers of particles. Thus, smearing is undone. In 

10 the concrete, the parameters are chosen such that a sum is minimized of the 
squares of the differences between measured numbers of reflections with 
amplitudes in different amplitude ranges and predicted numbers. Such a 
sum is an indication for the probability of the combination of measured 
reflections, but of course also other indications of the probability can be 

15 used, such as the product of the probabilities of the different measured 
amplitudes. 

A prerequisite for the technique used is that the reflections of 
individual particles can be distinguished. If reflections of several particles 
are measured in<hscriminately, a part of the reflections will be masked. 

20 Thus, reflection of a larger particle can make reflection of a smaller particle 
invisible, but reflection of a smaller particle can also mask reflection of a 
larger particle if the larger particle is sufficiently much farther from the 
center of the beam than is the smaller particle. 

Accordingly, the technique gives reliable results only when 

25 concentrations are sufficiently low. The maximum usable concentration can 
be raised by minimizing the volume in which particles are indistinguishably 
measured, for instance by the use of a focused beam and small time 
windows in which echoes are accepted. However, there are limits to the 
applicability of such techniques. That is why the requirement of separate 

30 observation of individual particles limits the applicability of the technique. 
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It is an object of the invention to make the use of ultrasonic 
measurements for the characterization of particle size distribution in liquids 
useful for higher concentrations. 
5 The present invention is based on the realization that it is possible 

with a statistical model to correct for the effect of particles being missed and 
that consequently the known measuring technique can also be used for 
larger concentrations. 

The invention provides a method according to claim 1. In the method, 

10 the parameters of the particle size distribution are estimated with a 

maximum likelihood technique, using an expression for the probability of 
measurement of reflection amplitudes in terms of two factors. The first 
factor expresses the probability Po(A) of a reflection measurement of which a 
reflection with the value A forms part, independently of the question 

15 whether this value is masked in that also a reflection with a greater value 
occurs. This first factor basically corresponds to the whole expression for the 
probability used in the state of the art. 

The second factor in the expression used according to the invention 
comprises the probability that there is not also a reflection with a 

20 dominating value of the property, which would mask the measured value, 
forming part of the reflection measurement. Thus, the effect of masking is 
taken into account in the estimation of the parameters. In this way, the 
maximum likelihood estimation takes into account both the effect of 
spreading and the effect of masking jointly . 

25 

These and other objects and advantageous aspects of the invention will 
be described with reference to the following figures: 

Figure 1 shows an apparatus for characterizing a liquid 
30 Figure 2 shows numbers of reflections with different amplitudes 
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Figure 3 shows a reconstructed particle size distribution 
Figure 4 shows a calculated concentration 

Figure 5 shows a flow chart of a method of estimating parameters of a 
particle size distribution. 

5 

Figure 1 shows an apparatus for characterizing a liquid. The apparatus 
comprises a liquid channel 11, a control unit 10, a signal generator 12, an 
ultrasonic transducer 14, a detector 16 and a calculating unit 18. Control 
unit 10 is coupled to signal generator 12, detector 16 and calculating unit 

10 18. Signal generator 12 has an output coupled to transducer 14. Detector 16 
has an input coupled to transducer 14. Detector 16 has an output coupled to 
calculating unit 18. Transducer 14 is coupled to channel 11 for generating 
an ultrasonic beam. 

In operation, a liquid with particles floating therein flows through 

15 channel 11. Control unit 10 triggers consecutive reflection measurements in 
each of which signal generator 12 generates a pulse which is converted by 
transducer 14 to an ultrasonic wave in channel 11. Transducer 14 captures 
reflections of the pulse and feeds resulting signals to detector 16. Detector 
16 selects signals which occur in a particular time window with respect to 

20 the generated pulse. As a result of the sound velocity in liquid, the time 
window defines a zone from where particles can give rise to reflections. 
Detector 16 measures the amplitude of the reflection signal (if present) in 
the time window and feeds the measured value of the amplitude to 
calculating unit 18. Calculating unit 18 collects histogram information of 

25 the amplitudes. To that end, the calculating unit utilizes a number of 
amplitude ranges and cumulates counts of the number of reflection 
measurements that have yielded reflection amplitudes in the respective 
ranges. After a large number of reflection measurements have thus been 
performed, calculating unit 18 calculates further information about the 

30 particle size distribution from the histogram information. 
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Fig. 2 shows a simulated example of histograms for liquids in which 
different concentrations of particle are present, while the particle size 
distribution fo(D) in each case is the same (that is, the fraction of the 
particles having a size between D and D+dD in each case is fo(D)dD for each 
5 value of D). It can be seen that, besides an increase of the numbers of 
reflections, also the shape of the amplitude distribution changes with 
increasing concentration. 

Fig. 3 shows, in dotted lines, estimates of the particle size distribution 
that can be calculated from the amplitude histograms on the basis of the 

10 known technique, for a number of different concentrations. Also shown is 
the true distribution. The known technique does not take into account the 
effect of the concentration on the shape of the amplitude distribution. The 
result is that the estimated particle size distributions are mutually different 
and at higher concentrations deviate from the true particle size distribution. 

15 The invention provides an improvement of the known technique of 

estimating the particle size distribution, so that the effect of the 
concentration of the particles can be corrected for as well. This improvement 
is based on the realization that the deviation of the shape of the histograms 
of Fig. 2 is a result of masking of reflection amplitudes which occur in the 

20 reflection measurement but are missed in the count because in the same 
reflection measurement one or more reflections with a greater amplitude 
occur. 

To that end, an expression is used for the probability distribution 
P(A)dA that a reflection measurement yields a measured amplitude which is 
25 in an infinitesimal interval between A and A+dA. This expression is in the 
form of 



P(A)=P 0 (A) Q (A) 
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Herein, P 0 (A) corresponds to the unmasked probability distribution, that is, 
probability of presence of a reflection, regardless of the question whether it 
is masked by another reflection having a greater amplitude. The factor Q(A) 
expresses the probability that in the same reflection measurement no 
5 particles cause reflections of an amplitude greater than or equal to A. 

The factor Po(A) corresponds to the probability distribution that is used 
in the state of the art for the whole probability distribution P(A) that an 
amplitude A in the interval A to A+dA is in fact measured (i.e. with neglect 
of the possibility that a reflection is masked). This factor P 0 (A) can be 
10 expressed in the particle size distribution fb(D) in a manner known per se. 
An expression for this factor is, for instance, 

Po (A)= C JdD fo(D) G(A i D) V m eas(D) 

15 Herein, C is the concentration of the particles (the average number of 
particles per unit volume), G(A | D) is the conditional probability that a 
reflection of amplitude A is detectable, if a particle of a size D yields a 
detectable reflection. The conditional probability G(A I D) expresses the 
effect that a particle of a size D gives rise to different reflection amplitudes 

20 at different distances from the center of the beam. V m eas(D) is the total 
volume of liquid from where a particle of size D can yield a detectable 
reflection. This volume is defined in depth by the time window from where 
reflections are accepted during the measurement and laterally by the decay 
of the beam intensity to a point where the amplitude of the generated 

25 reflections no longer projects above the noise level A n . 

In the frequently occurring case of a Gaussian distributed beam, the 
volume Vmeas(D) can be expressed as 

V meas (D)=7i: Az log(Ao(D)/ A n )/2K 



30 
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Herein, Az is the depth interval from where reflections are accepted, k is the 
standard deviation which is determined by the rate at which the intensity of 
the beam decays as a function of the distance to the center of the beam. 
Ao(D) is the maximum amplitude that a reflection of a particle of size D 
5 would yield if the particle were located in the center of the beam. 

For a Gaussian beam, also an expression for the conditional probability 
G(A | D) can be given: 



G(A|D)= l/{A*log(Ao(D)/A n )} 

10 

This expression applies provided that Ao(D)>A>A n ♦ For amplitudes greater 
than the maximum amplitude Ao(D) attainable with a particle size D, or 
smaller than the noise level A n , the conditional odds of detection G(A | D) 
equal zero. 

15 It will be clear, incidentally, that for beams having a distribution other 

than Gaussian, correspondingly different expressions for V m eas(D) and 

G(A | D) will be used. 

The newly introduced factor Q(A), which expresses the effect of 

reflections being masked by other reflections having an amplitude greater 
20 than or equal to the amplitude A, can also be expressed in terms of the 

concentration C and the particle size distribution fo(D). Normative for this 

is the integral expression 



W (A) = JdA'JdD &(D) G(A* | D) 

25 

Herein, the integral over the amplitude is taken from the noise level up to 
the amplitude A, for which the odds of detection are being calculated. The 
factor Q can be approximated by 



30 



Q(A)= exp[-C(l-W(A))] 
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Or, in approximation, by a factor 

Q (A) = [ 1-p +p (1-W(A)/Vmeas) ] n 

5 wherein the sum over k is taken from 1 to n; n is an integer not greater than 
the maximum number of particles that physically fit into the volume from 
where measurable reflections occur, but amply greater than the average 
number of particles that yield measurable reflections and p=CVmeas/n, 
where Vmeas is the volume from where the largest particles can still 
10 generate echoes above the noise level. 

An alternative expression for the P 0 (A) factor is for instance 

Po (A)= C JdD fo(D) 5V(A,D)/ dA 

15 Herein, C is the concentration of the particles (the average number of 

particles per unit volume). V(A,D) is the volume from where a particle of 
size D yields an echo of an amplitude greater than or equal to A. This 
volume is defined in depth by the time window from where reflections are 
accepted during the measurement and laterally by the decay of the beam 

20 intensity with the distance to the center of the beam. 

In the frequently occurring case of a Gaussian distributed beam, the 
volume V(A,D) can be expressed as 

V(A, D)=tc Az log(Ao(D)/ A)/2k 

25 

Herein, Az is the depth interval from where reflections are accepted, k is a 
parameter which characterizes the rate at which the intensity of the beam 
decays as a function of the distance to the center of the beam. Ao(D) is the 
maximum amplitude that a reflection of a particle of size D would yield if 
30 the particle were located in the center of the beam. It will be clear, 



WO 2005/057183 



PCT/NL2004/000857 



9 



incidentally, that for beams having a distribution other than Gaussian, 
correspondingly different expressions for V(A, D) will be used. 

The newly introduced factor Q(A) expresses the probability that 
besides the echo of amplitude A no echoes of greater amplitude occur. This 
5 describes the effect of reflections being masked by other reflections of an 
amplitude greater than or equal to the amplitude A. The factor Q(A) can 
also be expressed in terms of the concentration C and the particle size 
distribution fD(D). Normative for this is the integral expression 

10 W (A) = JdD &(D) V(A,D) 

The factor Q can be approximated by 

Q(A)= exp[-CW(A)] 

15 

On the basis of the expression P(A)= P 0 (A) Q for the probability distribution 
of measurements of reflections of amplitude A that are not masked by 
reflections of greater amplitude, the particle size distribution fb(D) is 
estimated from the system of reflection amplitudes that has been measured 
20 in a large number of measurements. This is done, for instance, by looking 
for an estimate fo(D) which minimizes the sum E of the squares of 
differences between the predicted numbers of measurements N pr ed of 
amplitudes in a series of amplitude intervals IAi and the actually measured 
numbers Nmeas 



25 



E = li (NmeasOAi)- N P red(IAi))2 



30 



This sum is in fact an indication for the probability of the measured system 
of amplitudes, so that the particle size distribution &(D) is chosen that 
makes this system most probable. Of course, instead of using the sum of the 
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squares, maximization can also be done using other indications of the 
probability, for instance using a product 



lb Po(Ai) Q (Ai) 

5 

of the probabilities of the different measurements Ai. 

In this estimation, preferably a parameterization of the particle size 
distribution fb(D) is used, for instance in the form of a multimodal 
distribution 

10 

fo(D) = li ci exp(-(D-Di)2/2di2 



(i=l, 2..M). In this case, the parameters ci Di and ai 2 are estimated such that 
according to the indication used, the probability of the measured system of 

15 reflection amplitudes is maximized. Techniques for such maximizations are 
known in general, and by applying them to an expression Po(A) Q (A) for the 
probability of measured (and not masked) reflection amplitudes, an estimate 
for the parameters is obtained. 

Figure 3 further shows different particle size distributions which are 

20 thus calculated from the numbers of reflections in the examples of Fig. 2. 
These distributions virtually coincide, within the accuracy of the figure, 
with the true distribution. It will be clear that thus a better estimate of the 
particle size distribution is obtained. 

Figure 4 shows the consequences of the estimate of the particle size 

25 distribution for the concentration. The concentration is the fraction of the 
liquid that is formed by the particles. Plotted horizontally is the true 
concentration and plotted vertically is the concentration calculated from the 
(simulated) measurements. The solid line indicates where calculated and 
true concentration are equal. The circles indicate the concentration 

30 calculated with the new technique and the crosses indicate the 
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concentration calculated with the state of the art. The figure shows that 
especially for higher concentrations clearly better results are obtained. 

It will be clear that the invention is not bound to this specific manner 
of estimating the particle size distribution. Thus, for instance, a different 
5 form of parameterizing the particle size distribution can be used, for 

instance on the basis of the values (and/or derivatives) of the particle size 
distribution for a number of particle sizes and an interpolation of the 
particle size distribution between these particle sizes, or, for instance, in the 
form of a quotient of two polynomials, whose coefficients form the 

10 parameters. Also, approximations for the different terms in the 
mathematical formula for the probability P(IA) can be used. 

Further, for instance, first an estimate can be made of the factor Q(A) 
as a function of the amplitude A (for instance by provisionally estimating 
the particle size distribution and subsequently calculating the factor with 

15 that provisionally estimated distribution), and then the measured numbers 
of reflections N m eas(A) with different amplitudes can be corrected with this 
factor Ncorr(A)= N m eas(A)/Q(A). Next, from the corrected numbers, the 
particle size distribution can be estimated in the manner known from the 
state of the art. 

20 Although the invention has been described in terms of the explicit 

particle size distribution fb(D), it will be clear that for certain applications it 
is not necessary to actually calculate this particle size distribution explicitly. 
When only a quantity such as for instance the mass density of particles is 
needed, which follows from the integral 

25 

JdD fo(D) m(D) 



then this quantity, which is a parameter of the particle size distribution, can 
also be calculated implicitly when estimating the parameters of the particle 
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size distribution that makes the system of measured reflection amplitudes 
the most probable. 

The technique used is not limited to measurements on amplitudes, but 
can be extended to any form of measurements of properties of the reflections 
5 whereby, in the presence of more particles causing reflections, the effect of a 
dominating particle masks the effect of the other particles. 

The technique is especially also applicable to liquids that contain 
different types of particles, when a distinction can be made as to which type 
of particle gives rise to a reflection measurement. This distinction may for 

10 instance be visible on the reflection signal itself, for instance by the sign of 
the first peak in a received ultrasonic reflection signal in response to an 
ultrasonic pulse of a particular sign. In the case of water in which oil 
droplets and sand grains are dissolved, for instance, reflections on sand 
grains and oil droplets can be distinguished by the sign of the reflection. 

15 The corrected probability distribution Pj(A)dA for a measurement of a 

reflection of a particle of type j, with an amplitude in a range between A and 
A+dA, can in this case be described with 



Pj(A)=P j0 (A) Ili Qi (A) 

20 

Herein, different factors Qi (A) describe for each different type of particles i 
the probability that no particles of the type i lead to detection with greater 
amplitude than A. The factors Q i (A) are each expressed in the above- 
described manner in the particle size distribution &i(D) of a respective type 
25 of particles i, and so are the probabilities Pj 0 (A). With these expressions for 
Pj(A) for different types of particles and a system of measurements of 
reflection amplitudes that are classified according to the type of particles 
involved, parameters of the particle size distributions fDi(D) can be 
estimated in the above-described manner. Thus, for instance, particle size 
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distributions can be estimated which minimize a composite sum of squares 
of deviations 

E = Zi Im (Ni m eas(IA m )- Ni pre d(IA m ))2 

5 

of the numbers Nimeas of measured reflections of particles of type m in 
different amplitude ranges IA m and the predicted numbers Ni pr ed of 
measured reflections in these ranges. 

Fig. 5 shows in summary the method for estimating the particle size 

10 distribution. In a first step 51, a pulse is generated and in a second step 52 
it is detected whether in a particular time window after the pulse a 
reflection above the noise level returns. In a third step 53 it is determined in 
which of a number of intervals IA the amplitude falls and in a memory a 
count for the interval in question is increased. From a fourth step 54, the 

15 previous steps are repeated until a particular number N of reflection 
experiments have been performed. Next, a fifth step 55 is performed in 
which parameters of an initial particle size distribution are chosen. With 
these, in a sixth step 56, using the formula P(A)=P 0 (A) Q (A), a prediction of 
the measured numbers is calculated and a difference with the measured 

20 numbers is determined. If the difference is sufficiently small, the method 
stops in a seventh step 57. If not, in an eighth step 58, the parameters are 
adjusted such that the difference will expectably decrease, and the method 
is repeated from the sixth step 56. Of course, estimation in the method is 
preferably carried out with a computer. 

25 It will be clear that the invention is not limited to the method 

described. For instance, in the sixth step 56, the difference does not need to 
be calculated explicitly, but instead information can be calculated that is 
needed to choose the adjustment of the parameters, and the method can 
stop in the seventh step 57 if the adjustment is less than a threshold. 

30 



